Tb 2 Ti 2 07, a two singlet spin-liquid pyrochlore 
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The pyrochlore compound Tb2Ti207 is known to remain in a spin-liquid state down to the 
lowest attainable temperature (0.05 K), whereas current theories predict it should order into an 
antiferromagnetic structure. A number of models have been developed in order to resolve this 
discrepancy, but they could not obtain a spin-liquid ground state. We present here an explanation 
for the lack of magnetic order in Tb2Ti2 07 based on the presence of a Jahn- Teller- like distortion 
from the local trigonal symmetry and on the physics of the two-singlet system coupled by exchange. 



Frustration in materials crystallizing in the pyrochlore 
lattice has been a subject of interest since the pioneer- 
ing work of J. Villain and an intense burst of activity 
in the field has taken place in the last ten years (for a 
review, see Ref.Q). Among the various compounds, the 
titanates R2Ti 2 07 (R=rare earth) occupy a special po- 
sition since many different ground states arc observed 
for the different R 3+ ions and an understanding of their 
basic properties is now achieved @). The Tb titanate 
Tb 2 Ti 2 07 stands as an exception, since its most salient 
property, i.e. the absence of magnetic ordering of the 
Tb 3+ moments remains unexplained. Pioneering neu- 
tron and muon measurements Q evidenced a cooper- 
ative paramagnetic state, where short range correlated 
Tb moments keep fluctuating down to at least 0.05 K, 
failing to develop long range order (LRO) in spite of a 
high Curie- Weiss temperature (~ —14 K) indicative of 
antiferromagnetic (AF) interactions. Since its discovery, 
Tb 2 Ti207 has remained a theoretical puzzle as well as 
an experimental challenge. It is considered as the best 
example of a spin liquid (SL) Q in the pyrochlore lat- 
tice, in which quantum fluctuations are at play through a 
peculiar crystal electric field (CEF) scheme of the Tb 3+ 
ion. Indeed, it shows a low lying energy level (a doublet 
at A ~ 18 K above the ground doublet) 0-0], in contrast 
with classical Ising spin ices with Ho or Dy 0] where 
A ~ 250 — 300 K. The paramagnetic spin correlations 
were understood taking into account these two low lying 
doublets Q, and low energy spin fluctuations were in- 
terpreted by a dynamically induced frustration involving 
virtual crystal field excitations, leading to the concept of 
" quantum spin ice" [To| • Indeed, local spin ice like config- 
urations have been observed under applied field [ill, [12 1 . 
In spite of all these attempts, the SL ground state of 
Tb 2 Ti207 remains unexplained so far. 

The Tb 3+ ion is a non-Kramers (or Van Vleck) ion, 
i.e. it possesses an even number (8) of electrons in its 
4/ shell. Its ground spin-orbit state has a total mo- 
mentum J=6 and a Lande factor gfj=3/2. A specific 
property of Tb2Ti 2 07 is that a strong 4/ electron-strain 
interaction combined with its peculiar CEF level struc- 



ture should lead to a cooperative Jahn- Teller (JT) transi- 
tion at Tjt ~0.1K (l3j |. which consists in a spontaneous 
symmetry lowering. Recently, high resolution x-ray scat- 
tering measurements down to 0.3 K [1J| have evidenced 
structural fluctuations below 15 K, which are thought to 
be precursor to such a Jahn- Teller transition and which 
would occur between cubic and tetragonal symmetry. 

We propose here that a cooperative tetragonal distor- 
tion develops in Tb 2 Ti 2 Cv at low temperature, which 
breaks the local trigonal symmetry {D^d) at the rare 
earth site. The degeneracy of the Tb 3+ CEF doublets 
is then lifted, and the low lying CEF states are two sin- 
glets separated by a small energy gap. For non-Kramers 
rare earths, the singlets are "non-magnetic", i.e. (J a )=0 
(a = x, y, z), in the absence of a magnetic/exchange field. 
Exchange is nevertheless present at low temperature and 
tends to establish a magnetic ground state (and thus a 
magnetic LRO) by mixing with excited states. This com- 
petition between crystal field and exchange can result in 
a paramagnetic ground state if the exchange is not strong 
enough |15| . In other words, there is a critical value of 
interionic exchange below which the material does not 
show magnetic LRO down to T=0. Short range and dy- 
namic spin correlations are present in this phase, which 
can be called a spin-liquid. We believe this mechanism is 
at play in Tb 2 Ti 2 07 and explains the absence of LRO. 
In this work, we apply this idea to explore in detail the 
consequences of a distortion from trigonal symmetry on 
the ground state of Tb 2 Ti207. This distortion probably 
occurs over a finite distance and during a finite lifetime at 
low temperature and preserves the overall cubic symme- 
try of the lattice. For zero or small distortion, we show 
that AF exchange would lead to a magnetic LRO ground 
state, but that above a threshold distortion value, there 
appears a range of exchange constants where no LRO 
can settle in, i.e. the system remains in a spin liquid 
state down to T—0. We also show that the assumption 
of isotropic exchange is not sufficient, and that the in- 
troduction of anisotropic exchange, known to be present 
in the pyrochlore titanates (lij . allows to quantitatively 
reproduce most properties of Tb2Ti2 07. 
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FIG. 1. (T,A) phase diagrams with the CEF level scheme 
of Tb2Ti2C>7, with isotropic exchange and dipolar interac- 
tion. Upper left panel: no distortion; lower panel: with 
a tetragonal distortion Dq =0. 25 K; upper right panel: SL 
phase limits as a function of the distortion stength Dq. The 
transition line to the AF phase is in green, to the OSI phase 
in red, and the SL phase is in blue. 



We start by computing the phase diagram as a func- 
tion of exchange. We consider the total hamiltonian for 
a Tb 3+ ion: T~L — Hcef + %ex + 'Hdip, containing the 
CEF interaction as in Rcf.[7], the short range exchange 
interaction with the 6 nearest neighbours and the infinite 
range dipolar interaction calculated using Ewald summa- 
tion method [17[. The ground state at a given tempera- 
ture is obtained by a self-consistent mean field calculation 
involving the 4 sites on a tetrahedron. Then, only mag- 
netic structures having k=0 propagation vector can be 
reached. The exchange interaction is described in terms 



of a molecular field tensor A, with J = ^ Mb , 

where J is the usual exchange integral, and the molecular 
field acting on a given ion is: H mo i = A J2k=i 6 m k- We 
first set zero distortion (upper left panel in FigJT]). Mag- 
netic ordering is found to occur for all A values. We find 
there is a crossover value: \ cr = — 0.10 T//zb below which 
the ground state is an AF state ("all in/ all out"), and 
above which it is an ordered spin-ice phase (OSI, "two in - 
two out") akin to that observed in Tb2Sn207 [1S|] . In the 
OSI phase, the Tb moments are slightly tilted with re- 
spect to their local [111] axis. On decreasing |A| towards 
the ferromagnetic region, they tend towards [00 11, recov- 
ering results similar to those obtained in Ref.[19|. As 
in the spin ice phase diagram of Ref . [2C| , we obtain the 
OSI phase for not too large antiferromagnetic exchange 
constants. 

We next introduce a tetragonal distortion Hq = 
DqJ^, where OZ is the cubic [001] axis. In the local 
frame with [111] as z-axis, the distortion term writes: 



U Q = ^ [2Jl + Jl + V2(J X J Z + J Z J X )\. (1) 
An interesting behaviour appears in the phase diagram 



as one switches on the distortion: there opens a gap in 
the low |A| region of the OSI transition line for Dq > 
0.22 K. The lower panel of Fig JT] shows this behaviour 
for Dq=0.25K, which is the strength of the distortion in 
Tb2Ti207, as will be shown below. For A values inside 
the gap, the spontaneous moment at each Tb 3+ site is 
zero (down to 0.05 K). This region corresponds to a phase 
with no long range order down to very low temperature 
and is analogous to that found in Ref. [3]. It should be 
the basic reason for the absence of LRO in Tb2Ti20y. 
The two OSI phases obtained on either side of the gap 
are slightly different: phase OSI I is similar to the OSI 
phase obtained with zero distortion, whereas phase OSI 
II has a different symmetry. The gap width increases as 
Dq increases above 0.22 K, as illustrated in the upper 
right panel of FigJTJ where the lower and higher limits of 
the SL phase |A_| and |A + | are plotted versus Dq. 

The distortion induces important changes in the eigen- 
functions of the ground states. In the basis |m) = | J = 
6; J z = m), the ground trigonal wave- functions are: 

~ a| - 4) + 6|5> and |^> -6| - 5) + a\4) (2) 

with a ~0.96 and b ~0.23. This non-Kramers doublet 
is such that the transverse matrix element (V^l^alV'i) i s 
zero in zero magnetic field and remains close to zero as 
a moderate field is applied. Introducing the distortion 
§T§ yields new eigenstates which are entangled singlets of 
symmetric/antisymmetric type: 
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A particular property arises in Tb2Ti 2 07 from the struc- 
ture of states ((2]): the tetragonal distortion lifts their 
degeneracy at first perturbation order, giving an energy 
separation of the two ground singlets 8 = 12abDQ\f\\. 
For D Q =0.25K, this yields 5=2.2 K=0.046THz. This 
property also holds in Tb2Sn207, but not in Ho2Ti207, 
for instance, where the splitting of the singlets for 
-Dq=0.25K would be a few mK. By contrast with 
the trigonal states ©, the transverse matrix element 
(t^s | |"0a) is non- vanishing and large. 

The low temperature inelastic neutron spectra, both 
in zero field @, and in finite field 21 , 22| , show a low 
energy inelastic line which gives a direct evidence for the 
presence of the distortion and allows its strength to be 
determined. At 1.6 K in zero field, this line is not well 
resolved from the quasi-elastic line, but it is visible at a 
low energy S ~ 0.04 THz besides the CEF transition near 
0.4 THz (see Figl2]a). Its intensity decays more rapidly at 
large Q than the square of the form factor, demonstrating 
its magnetic origin. At 0.4 K in an applied field of 3 T, the 
quasi-elastic line is much narrower and the inelastic line, 
which has slightly shifted to an energy 6 ~0. 055 THz, is 
very well resolved (see FigJ5]b). This line arises from a 
degeneracy lifting of the ground doublet, and the simu- 
lations in the lower panel of FigJ5] demonstrate that the 
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FIG. 2. Neutron inelastic spectra in Tb2Ti207. The upper 
panels show the data, with dashed Lorentzian-shaped lines 
marking the peak positions. The lower panels show the sim- 
ulations with the full hamiltonian. a: spectra in zero field 
at T—1.6K; data from Ref.Q and simulations with a tetrag- 
onal distortion Dq = 0.25 K (blue solid line) or with a Zee- 
man interaction with i/ cx =lT (red dashed line), b: in-field 
spectra with H=3T applied along [110], at T=0.4K; single 
crystal data from Ref . [22l[ and simulations with the distortion 
-Dq=0.25 K (blue solid line) or without distortion (red dashed 
line). 



tetragonal distortion alone can generate this degeneracy 
lifting and account for a sizeable intensity transition (blue 
solid line) . Both the hypothetical presence of a static ex- 
change field of 1 T for the zero field spectrum, or the 
presence of the applied field alone for the 3 T spectrum 
result in a quasi-extinct line (red dashed lines in FigJ5]). 
These extinction rules originate in the above consider- 
ations about the wave-functions, the intensity iy of a 
transition between states \i) and \j) being proportional 
t° Eb=i s z K*l^a|i)| 2 < From the position of this line, 
the strength of the tetragonal distortion in Tb 2 Ti207 is 
estimated to be Dq=0.25K. 

The low energy part of the CEF spectrum in the pres- 
ence of the distortion consists in 4 non-degenerate states, 
i.e. 2 sets of closely spaced singlets with a splitting of 
about 2K, separated by 18 K. Therefore, specific heat 
measurements carried out from a temperature lower than 
the splitting of the ground singlets should derive an en- 
tropy release AS close to i?ln4 as soon as the 4 states 
are populated, instead of R In 2 for 2 degenerate doublets. 
Specific heat data in Tb2Ti2C>7 23j, |24j have indeed ob- 
tained AS ~ i?ln4 between 0.15 and 20 K, confirming 
the level structure proposed by our model. 

One must now examine whether the other available 
experimental data in TL>2Ti207 are compatible with A 
values within the SL gap associated with Dq=0.25K, 
i.e. -0.028 < A < -0.015 T/fi B - From the high tem- 
perature susceptibility data, an estimation was given: 
A ~ -0.05T///B 0, not far from, but outside the SL 
range. However, other data (the thermal variation of the 
local susceptibility [l6| and the field evolution of the in- 



FIG. 3. (A c , A a ) phase diagram with the CEF level scheme 
of Tb2Ti207 at T=0 in the presence of anisotropic exchange, 
dipolar interaction and a tetragonal distortion Dq=0.25K, 
with the fixed value Af, = — 0.13T//iB- The red region 
has the OSI phase as ground state, the blue region the 
SL phase and the green region the AF phase. The point 
Xo ={-0.045, -0.13, -0.07}T/|LtB, marked with a pink star, 
could correspond to Tb2Ti2C>7 (see text). 
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FIG. 4. Calculations of the field variations, for H // [110], of 
the a and /3 Tb 3+ moments at 1.6 K, showing the spin-melting 
field for f3 moments, obtained with the molecular field tensor 
Ao (see Fig[3]) (solid lines). The dashed lines represent the 
calculation with the isotropic molecular field constant A = 
-0.025 T//xb. The data are from Ref. [H. 



duced magnetic structure [12|) point to the presence of 
anisotropic exchange. So we have determined the phase 
diagram at T—0 assuming a 3-component molecular field 
tensor {A a ,Af,,A c }, where the c axis is chosen along the 
bond linking two Tb ions. A representative cut of this 
phase diagram, for A^ = — 0.13T//ib, is shown in FigJ3] 
For large |A| values, the OSI phase is the ground state 
(in red) , but a sizeable SL region (in blue) together with 
an AF region (in green) are found for lower |A| values. 

A quantity which is rather sensitive to the A values 
is the in- field behaviour of the magnetic structure flij ]. 
In particular, the /3-moments (those having their trigo- 
nal axis perpendicular to the field H applied along [110]) 
undergo a "spin-melting" at H m = 1.5(5) T, where their 
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FIG. 5. Thermal variation of the single crystal magnetic sus- 
ceptibility calculated with the CEF parameters for Tb2Ti207, 
with the anisotropic Ao tensor derived above, in the absence 
and in the presence of the tetragonal distortion. Insert: in- 
verse powder susceptibilities with the same assumptions and 
legend, and experimental data from Ref. 25] (blue squares). 



magnitude vanishes, as shown in FigSl The curves com- 
puted using the anisotropic Ao tensor in the SL region of 
FigEl (solid lines) reproduce the data better than those 
obtained using the isotropic value A = — 0.025 T//ljb 
(dashed lines), corresponding to a point in the SL range 

of Figru 

The magnetic susceptibility x(T) at very low temper- 
ature should reflect the presence of the distortion, as 
show the calculated curves in Fig [5] for Dq—0 (black 
curves), x(T) shows a Curie- Weiss behaviour, whilst for 
Dq=0.25K (red curves) it shows a Van Vleck-like be- 
haviour with a saturation below IK, as expected for a 
singlet ground state separated from the excited state by 
about 2K. The powder susceptibility measured by mag- 
netometry 25] does not saturate (see insert of Figj5l blue 
squares), but l/\ computed without distortion is farther 
from experimental data than that computed with distor- 
tion. The powder susceptibility can also be measured 
through the /iSR Knight shift in a transverse field. The 
Knight-shift data of Ref. 1261 d o not evidence any net sat- 
uration, but those in Ref.[27| do show a clear saturation 
below 1-2 K, in good agreement with our model. 

The numerous neutron scattering studies performed in 
Tb 2 Ti 2 7 (EH El have shown that the Tb 3+ ions be- 
have like short range correlated moments down to low 
temperature. Our picture is fully compatible with fluc- 
tuating Tb moment components with sizeable values. 
The fast spin fluctuation rate down to the 0.05 K range 
derived by ^SR and neutron spin echo [28] can be 
explained by exchange/dipole driven relaxation between 
Tb 3+ singlets of type ([3]), which are linked by large ma- 
trix elements of J z . Indeed, according to the Fermi rule, 
the relaxation rate is proportional to \(ip a \ J z \^s}\ 2 - Bc- 



cubic cells, with the same propagation vector k=(001) as 
the LRO induced under stress 13111. 



As a conclusion, we think the absence of magnetic LRO 
in Tb2Ti207 is caused by the development at low temper- 
ature of a precursor Jahn- Teller distortion which makes 
this material a two-singlet system, where AF exchange is 
not strong enough to allow for magnetic ordering. Tak- 
ing also into account exchange anisotropy, most of the 
properties in Tb2Ti207 can be explained. The in-field 
spin-wave spectrum in Tb2Ti207 [2l[ will be accounted 
for in a future work. Our model can also be applied suc- 
cessfully to explain the characteristics of the OSI phase 
in Tb2Sn 2 07 [lSj] ; this will be the subject of a future 
publication. 



low about 0.5 K, the spin fluctuations freeze at the 10 s 
time scale [30j j , and coexist with spin-glass irreversibilities 
[25| and with magnetic correlations extending over a few 
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